Optimized optomechanical crystal cavity with acoustic radiation shield 
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We present the design of an optomechanical crystal nanobeam cavity that combines finite-element 
simulation with numerical optimization, and considers the optomechanical coupling arising from 
both moving dielectric boundaries and the photo-elastic effect. Applying this methodology results 
in a nanobeam with an experimentally realized intrinsic optical Q- factor of 1.2 x 10 6 , a mechanical 
frequency of 5.1 GHz, a mechanical Q-factor of 6.8 x 10 5 (at T = 10 K), and a zero-point-motion 
optomechanical coupling rate of g = 1.1 MHz. 



The use of radiation pressure forces to control and 
measure the mechanical motion of engineered micro- and 
nanomechanical objects has recently drawn significant 
attention in fields as diverse as photonics pQ, precision 
measurement [2], and quantum information science [3]. 
A milestone of sorts in cavity circuit- and optomechan- 
ics is the recent [4j [5] cooling of a mechanical resonator 
to a phonon occupancy (n) < 1 using cavity-assisted ra- 
diation pressure backaction [6j [7j. Backaction cooling 
involves the use of an electromagnetic cavity with reso- 
nance frequency uo Q sensitive to the mechanical displace- 
ment, x, of the mechanical resonator. The canonical sys- 
tem is a Fabry-Perot cavity of length L with one end 
mirror fixed and with the other end mirror of mass m e ff 
mounted on a spring with resonance frequency uj m . The 
coupling between the electromagnetic field and mechan- 
ics is quantified by the frequency shift imparted by the 
zero-point motion of the mechanical resonator, given by 
g g OM y / h/2m e ^uj rn , where g OM du Q /dx uj /L. 

One of many technologies recently developed to make 
use of radiation pressure effects are optomechanical crys- 
tals (OMCs) [8j [9]. Optomechanical crystals, in their 
most general form, are quasi-periodic nanostructures in 
which the propagation and coupling of optical and acous- 
tic waves can be engineered. In this work we present the 
comprehensive design, fabrication, and characterization 
of a quasi- ID OMC cavity formed from the silicon de- 
vice layer of a silicon-on-insulator (SOI) microchip. Our 
design incorporates both moving-boundary and photo- 
elastic (electrostriction) radiation pressure contributions, 
and simultaneously optimizes for optical and acoustic pa- 
rameters. 

The nominal unit cell of a nanobeam OMC, geomet- 
rically a silicon block with an oval hole in it, is shown 
schematically in Fig. [IJl The corresponding optical and 
mechanical bandstructure diagrams are shown in Figs.[]J 
and d, respectively. As indicated by the gray shaded re- 
gion in the photonic bandstructure, the continuum of un- 
guided optical modes above the light line precludes the 
existence of a complete photonic band gap (only a quasi- 
bandgap exists for the guided modes of the beam). The 
physical dimension of the unit cell block (see caption of 
Fig. Ilk) are chosen to yield a photonic quasi-bandgap sur- 



rounding a wavelength A ~ 1550 nm. The corresponding 
mechanical bandstructure has a series of acoustic bands 
in the GHz frequency range (the ratio of the optical fre- 
quencies to that of the mechanical frequencies is roughly 
the ratio of the speed of light to sound in silicon). By 
classifying the acoustic bands by their vector symmetry, 
we can again define a quasi-bandgap for the mechanical 
system in terms of modes of common symmetry (modes 
of different symmetry may couple weakly due to sym- 
metry breaking introduced by the fabrication process, an 
issue we discuss and mitigate later in the design process) . 

Using previously developed design intuition [9], we 
choose to focus on the optical "dielectric" band at the 
X-point and the mechanical "breathing mode" band at 
the T-point (these bands emphasized by thicker lines in 
Figs.[lJ and e). Fig{l]i and f shows the optical dielectric 
band X-point and the breathing mode T-point frequency 
as the unit cell is transformed smoothly from that in 
Fig.[l^ to that in Fig.[l]3. We see that such a transition — 
simultaneously reducing the lattice constant a and de- 
creasing the hole aspect ratio, = h y /h x — causes an 
increase in the X-point optical frequency and a reduction 
in the T-point mechanical frequency, pushing both opti- 
cal and mechanical modes into the quasi-bandgap of their 
respective bands (shown as shaded regions in Figs.JT^ and 
e). 

An optomechanical cavity can be formed by transi- 
tioning from the nominal unit cell in Fig. [TJi, forming the 
"mirror" region, to the defect unit cell of Fig. ^p. This 
cavity can be parameterized by the maximum change in 
a (denned as d = max{l — a/a n0 minai}) 7 the number of 
holes in the defect region, the maximum curvature of 
^/^nominai (plotted in Fig|2^i) and the of the center 
hole. Including the unit cell geometric parameters a, w, 
h x , and h y (t is fixed to 220 nm), an entire nanobeam de- 
sign is specified by these 8 values. Finite-element-method 
(FEM) simulations of a complete structure are used to 
determine the fundamental cavity mode frequencies (uo 
and o; m ), motional mass m e ff [10], and radiation-limited 
optical Q-factor, Q Q . To quantify the coupling rate, we 
consider both the frequency shift due to the moving di- 
electric boundary [11 and the photo-elastic effect [12], 
so that g OM = #om,mb +#om,pe- The moving boundary 
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FIG. 1: The (a) nominal unit cell with (a,t,w,h x ,h y ) = 
(436,220,529,165,366) nm and (b) defect unit cell with 
(a,t,w,h x ,h v ) = (327,220,529,199,170) nm of the OMC 
nanobeam cavity. The (c) optical and (e) mechanical band 
structure for propagation along the x-axis in the nominal unit 
cell, with quasi-bandgaps (red regions) and cavity mode fre- 
quencies (black dashed) indicated. In (c), the light line (green 
curve) divides the diagram into two regions: the gray shaded 
region above representing a continuum of radiation and leaky 
modes, and the white region below containing guided modes 
with y-symmetric (red bands) and y- antisymmetric (blue 
bands) vector symmetries. In (e), modes that are y- and 
z-symmetric (red bands), and modes of other vector symme- 
tries (blue bands) are indicated. The bands from which the 
localized cavity modes are formed are shown as thicker curves. 
Tuning of the (d) X-point optical and (f) T-point mechani- 
cal modes of interest as the unit cell is smoothly transformed 
from the nominal to the defect unit cell. 



contribution is given by [TP] , 

u £(Q-n)(AeEf, -Ae-^l) dS 



#OM,MB 



(1) 



where Q is the normalized displacement field 
(max{|Q|} = 1), n is the outward facing surface normal, 
E is the electric field, D is the displacement field, 
the subscripts || and _L indicate the field components 
parallel and perpendicular to the surface, respectively, 
e is the material permittivity, Ae = ^silicon — £ain an d 
^air ID- A similar result can be de- 
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FIG. 2: (a) Plot of the unit cell lattice constant, a, along the 
length of the nanobeam. (b) The normalized optical E y field 
and (c) the normalized mechanical displacement field |Q| of 
the localized optical and mechanical modes, respectively, (d) 
The normalized surface density of the integrand in eq. (IT]), 
showing the contributions to ^om,mb- (e) The normalized 
volumetric density of the integrand in eq. (|3|, showing the 
contributions to goM,PE- (f) Scanning electron microscope 
(SEM) image of the experimentally realized cavity. 



perturbation theory, 



#OM,PE 



u Q (E\£\E) 
2 /E-DdV" 



(2) 



where a is a generalized coordinate parameterizing the 
amplitude of Q. In an isotropic medium with refractive 
index n, we have = —Son 4 pijkiSki where p is the 
rank-four photo-elastic tensor and S is the strain ten- 
sor [13] . For silicon, a cubic crystal with point symme- 
try group m3m, with the x-axis and 7/-axis respectively 
aligned to the [100] and [010] crystal direction, we can use 
the contracted index notation with some simplification to 
write 



- eon 4 J 2Re{E*E y }p i4 S i 

+ 2Re{E*E z }p 44 S 5 + Re{E;E z }p M S 6 
+ \E x \ 2 ( Pu S 1 + Pl2 (S 2 + S 3 )) 
+ \E y \ 2 ( Pll S 2 +p 12 (S 1 + S 3 )) 

+ |-Ez| 2 (Pll^3+Pl2(^l+^2))l^ (3) 



Ae" 1 

rived for the photo-elastic contribution from first-order where (pu,pi2,P44) = (—0.094,0.0 1 7,-0.051) [T2l IT3]. 
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FIG. 3: (a) The unit cell of the phononic shield with param- 
eters (c a ,c h ,ct,t) = (534,454,134,220) nm. (b) Full in-plane 
mechanical band diagram with complete phononic bandgap 
(red region) and frequency of the acoustic mode of the cavity 
indicated (black dashed line), (c) SEM image showing the 
phononic shield (green) around the nanobeam. 



In order to optimize the nanobeam OMC design 
we have chosen to assign a fitness value F = 
— g • min{Q , Qcutoff }/Qcutoff to the simulations. The ad- 
ditional Qcutoff term (set to 3 x 10 6 ) is used to avoid 
unrealizably high simulated radiation-limited Q Q values 
from unfairly weighting the fitness. With this choice the 
nanobeam design has been reduced to an 8 parameter 
optimization problem amenable to a variety of numeri- 
cal minimization techniques. For a computationally ex- 
pensive fitness function with a large parameter space, a 
good choice of optimization algorithm is the Nelder-Mead 
method [14] ; as a downhill simplex method (compared to 
a gradient descent method) the search technique is resis- 
tant to simulation noise and discontinuities. To mitigate 
the problem of converging on local minima, the optimiza- 
tion procedure is applied to many randomly generated 
initial conditions. The resulting optimized nanobeam 
OMC design is shown in Fig. [2j with a simulated uj /2tt 
of 194 THz, radiation-limited Q Q of 2.2 x 10 7 , o; m /27r of 
5.7 GHz, and ra e ff of 127 fg. The total optomechanical 
coupling rate, g/27r, is 770 kHz, composed of a —90 kHz 
contribution from the moving boundary and a 880 kHz 
contribution from the photo-elastic effect. 

Unaddressed so far in our design is mechanical losses, 
a source of which is coupling (made unavoidable by 
inevitable symmetry breaking during the fabrication 
process) to unconfined modes of alternate symmetries 
present in the quasi-bandgap of the patterned nanobeam 
(Fig. [IJ). These acoustic radiation losses can be signif- 
icantly suppressed by surrounding the entire nanobeam 
inside a second acoustic shield consisting of a pattern- 
ing with a complete phononic band gap. The 2D cross 
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FIG. 4: (a) Normalized optical transmission spectrum, cen- 
tered at 1544.8 nm, showing the fundamental optical cav- 
ity mode of the nanobeam, with a measured intrinsic Q ,i — 
1 .22 x 10 6 . (b) Optically transduced thermal noise power spec- 
tral density centered at the mechanical frequency, o; m /27r = 
5.1 GHz, of the breathing mode, taken at T s — 10 K with 
the input laser red-detuned (A = cj m ; red curve) and blue- 
detuned (A = — cj m ; blue curve) from the cavity, (c) Mea- 
sured cooperativity, C, as function of intracavity photon num- 
ber, ricav, for red detuning A = uj m . 



phononic crystal design [TO] has previously demonstrated 
this to great effect [15], and the wide tuneability of the 
acoustic gap essentially allows the mechanical loss en- 
gineering to be decoupled from the design of the co- 
localized cavity modes of the nanobeam. By adjusting 
c a , Ch and c t in the 2D unit cell (Fig. |3^i), it can be seen 
from Figj3]o that the complete bandgap can be tailored 
to be centered on uj m . 

The optimized nanobeam design was fabricated from 
the [001] silicon device layer of a silicon-on-insulator 
wafer from SOITEC (resistivity 4—20 Q • cm, device layer 
thickness 220 nm, buried-oxide layer thickness 2—3 |Ltm). 
The cavity geometry was denned by electron beam lithog- 
raphy followed by inductively-coupled-plasma reactive- 
ion etching to transfer the pattern through the 220 nm 
silicon device layer. The cavities were then undercut us- 
ing a 1:1 HF:H20 solution to remove the buried oxide 
layer, and cleaned using a piranha/HF cycle. The sili- 
con surface was hydrogen-terminated with a weak 1:20 
HF:H20 solution for chemical passivation [T6] . 

Characterization of the OMC nanobeam devices was 
performed in a continuous- flow helium cryostat, under 
vacuum and at a sample mount temperature of T s ^6K. 
A tapered optical fiber [17], positioned in the optical 
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near-field (~100 nm) using a set of low- temperature- 
compatible piezoelectric stages, is used to evanescently 
couple laser light into and out of the devices. A tun- 
able external cavity diode laser (New Focus, model 6728) 
is used to scan across the wavelength band from A = 
1520—1570 nm. The resulting normalized optical trans- 
mission scan of an OMC nanobeam cavity with resonance 
wavelength at A G = 1544.8 nm (uj /2tt = 194 THz) is 
shown in Fig. [4^i. A Lorentzian fit to the measured op- 
tical linewidth (5\ = 1.7 pm) yields a fiber-taper- loaded 
optical Q-factor of Q Q = 9.06 x 10 5 , which for the mea- 
sured on-resonance transmission of 55% corresponds to 
an intrinsic Q-factor of Qo,i = 1.22 x 10 6 . The corre- 
sponding total cavity (energy) decay rate, (bi-directional) 
fiber taper waveguide coupling rate, and intrinsic de- 
cay rate are k/2tt = 214 MHz, K e j2i\ = 55 MHz, and 
Kijln — 159 MHz, respectively. 

Spectroscopy of the mechanical mode is performed by 
tuning the frequency of the input laser {uj\) to either a 
mechanical frequency red- or blue-detuned from the op- 
tical cavity resonance (A = (uo — uji) = ±u; m ). The 
optomechanical backaction under such conditions results 
in an optically induced damping of 7om = ±4G 2 /ft for 
A = ±uj m [18] , where G = ^n CdiV g is the parametrically- 
enhanced optomechanical coupling. The intracavity pho- 
ton number can be related to the laser input power (Pi) 
through the relation, n cav = 

The thermal Brownian motion of the mechanical mode 
is imprinted on the transmitted laser light as a side- 
band of the input laser resonant with the optical cav- 
ity resonance. As shown in Fig. [4}d, the resulting in- 
tensity modulation of the photodetected signal gives rise 
to a Lorentzian response in the photocurrent electronic 
power spectrum around the 5.1 GHz resonance frequency 
of the breathing mechanical mode. The intrinsic me- 
chanical damping of the breathing mode is found by av- 
eraging the measured mechanical linewidths under red 
and blue detuning, 7^ = (7+ + 7-)/2 = 7.5 kHz, where 
7± = 7^ ± 7om- The corresponding intrinsic mechanical 
Q-factor is Q m ,i = 6.8 x 10 5 . Experiments with different 
surface preparations (with and without HF-dip prior to 
testing) indicate that Q m ^ is extremely sensitive to the 
surface quality, and is likely not limited by bulk mate- 
rial damping at these temperatures [15]. A plot of the 
cooperativity, C = 7oM/7i = 7+/7i — 1? is shown in 
Fig. [4)3 versus n c , yielding the zero-point optomechani- 
cal coupling rate of g = 1.1 MHz for the breathing mode, 
very close to the simulated value (discrepancy here is at- 
tributed to the uncertainty in the silicon photo-elastic 
coefficients, which were extrapolated from data at wave- 
lengths of 3.4 urn and 1.15 [im \VH\19\ ). 

These measured device parameters place the system 
in the deeply resolved sideband regime (uj m /K > 30), 
suitable for efficient laser cooling of the mechanical res- 
onator into its quantum ground state. Calibration of the 
thermal Brownian motion [5] at low incident laser power 



(^cav ^ 1) indicates that the local bath temperature of 
the 5.1 GHz breathing mode is X5 w 10 K, correspond- 
ing to a thermal bath occupancy of only n& = 43, and a 
thermal decoherence time of r t h = ^Qm/te^ ~ 0.5 M-s 
(>10 4 cycles of the mechanical resonator). With a mea- 
sured (simulated) granularity parameter of gj k% « 0.0069 
(0.12), such devices with further improvement in the opti- 
cal Q-factor, represent a promising platform for realizing 
quantum nonlinear optics and mechanics [20] . 
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